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PeoTRACL 
The objective of ee ech in the field of axially compressed 
cylinders is the establishment of sound design principles. Only 
through additional analytical and experimental effort can this goal 
be met. The most critical need at the present is believed to be 
for additional tests, which, unlike most previous tests, will be 
designed specifically to provide general information about the 
buckling process. This investigation adds to the general study on 
the buckling of thin-walled cylindrical shells fabricated from the 
polyester film, Mylar. It was concluded that the number of cir- 
cumferential buckles, which controls the postbuckling behavior of 
cylinders, is dependent on the geometric parameters, R/t and 
L/R. Such dependence is approximated by de Neufville's solution: 
we o.e [RE] 

; L/R 
It is also concluded, for the range of the R/t and L/R para- 
meters considered in this investigation, that the ratio of the buck- 
ling stress is: (1) independent of the parameter L/R, and (2) 
dependent on the parameter R/t, with increasing values of the 


ratio R/t causing a decrease in 0/o el” 
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mist OF SYMBOLS 
buckling load 
classical buckling load 
number of circumferential buckles 
radius of cylindrical shell 
effective length of cylinder 
wall thickness of cylinder 


Young's Modulus of material 


1. Introduction 

For over thirty years the serious discrepancy between 
theoretical and experimental results for the buckling of axially 
compressed cylindrical shells has been recognized, but the 
phenomemon has had little understanding. Figures 1 and 2 
are block diagrams depicting the chronology on the investiga- 
tion of buckling stress and buckling mode of thin cylindrical 
shells. In 1932 Fliigge calculated the theoretical buckling 


stress for asymmetric buckling as 0 , = 0. 6 Et/R. Experi- 


] 
mental investigations found that the ratio of the actual buckling 
stress to the classical buckling stress was quite low with great 
scatter in the test results. The study of Karman and Tsien in 
1941 presented a possible explanation of the scatter in the test 
results and the reduction in the buckling load [6]. By the use 

of the theory of finite deflections, they managed to show that 

the shell stability is quite weak at load levels considerably be- 
low the classical buckling load. The essential result of their 
investigation was the statement that the axially compressed 
cylinder, after having reached the critical buckling load must 
snap through into another state of equilibrium which is con- 
nected with a considerably smaller axial load [9]. Their theory 


was slowly refined in the following years. However, due to the 


assumption of an initially perfect cylinder, which was used in 


these studies, the highest load carried according to the non- 
eae theory still 1s the buckling load predicted by classical 
theonmy In 1950 eemnell and Wan suggested that the cause for 
reteentte : and reduction of buckling loads was due to initial 
imperfections in the circular cylinder [4]. Their investigation 
ar CWE that the actual buckling load is considerably reduced, 
compared with the classical buckling load, due to initial imper- 
fections in the order of only a few tenths of the thickness of the 
shell wall, and that the difference between the actual and class- 
ical buckling loads depends on the amplitude of these initial 
imperfections [9]. However, the actual shape of the imperfec- 
tions in their analysis is slneewas and the exact degree of the 
effect for various types of imperfeceaie is still to be determined. 
Results in recent years by Almoth, Holmes and Brush [1], 
Babcock and Sechler [2], and Tennyson [8] show that very high 
values of the buckling load can be obtained when initial irregu- 
larities are carefully minimized. This indicates strongly that 
when cylindrical shells are sufficiently free of imperfections 
and ae loaded, the classical buckling load can be closely 
obtained and a lower buckling load is not always inevitable. 
Whether a useful theoretical solution will ever be feasible 
remains to be seen, but even if cylindrical shells are designed 
by use of empirical methods, it is important that we understand 


what happens when the shell buckles. 
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The concept of non-linear buckling showed that a buckled 


state could be maintained at o €o0 P but failed to show how the 


Y 
cylinder buckled without passing through the condition that 

i The mechanism of the transition is now the main 
problem in shell buckling. 

Calculations which are presently made on the number of 
buckles formed around the circumference of an axially com- 
pressed cylinder are in variance with those observed experi- 
mentally. From published experimental data, the variation of 
the number of circumferential buckling lobes with the cylinder 
parameters R/t and L/R have been obtained. Hoff suggested 
that the number of buckling lobes are interrelated with the crit- 
ical compressive loading [5]. The nature of this dependency 
poses a new problem in the understanding of the mechanism of 
the buckling of axially loaded circular cylinders. 

It is important to be able to determine the number of buck- 
ling waves at failure, N, because it is an index of the buckling 
mode, which in turn controls the behavior and strength of the 
cylinder after buckling. The number of circumferential buckles 


can be determined from the linear Flugge equations and is 


approximately predicted by de Neufville's solution in 1965 


as [7]: 1/4 
n= 2.2 [BA] . 


De Neufville has done extensive investigations on Mylar 
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cylinders in the region for R/t greater than 800 and less than 
1600. The major aim of the present paper is to investigate and 
validate de Neufville's analytical formulation for values of the 
parameter R/t less than 800. 

The primary factors which have been considered in this re- 
port are the influence of L/R, R/t, and o/o 4 on the number 
of circumferential buckles. In this investigation over 300 tests 
were conducted on over 50 cylinders fabricated from the 
polyester film, Mylar. 

The first part of the test program was designed to deter- 
mine if a normal distribution of results would be obtained for 
nominally identical test models. 

The remaining part of the test program was designed to 
determine the effect of varying the geometric parameters, 

L/R and R/t, onthe number of circumferential buckles. 

It is hoped that this study will stimulate the interest of 
others to perform additional work needed to complement and 
extend the range of geometric parameters considered, as the 
under stants of the buckling phenomena is still far from com- 
plete. It is also hoped that the net result of these investigations 
will bring some order to the study of cylindrical shells in axial 


compression. 
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2. Model Material 

All of the cylindrical models used in the experimental pro- 
gram were fabricated from Type A Dupont Mylar Polyester Film. 
Mylar is readily available in thin sheets, of 0.00015" (.15 mil) 
to 0. 014'' (14 mils) in thickness, and it is inexpensive. Extra- 
ordinary folding endurance is a characteristic over the full range 
of temperature and humidities. Mylar was used because of its 
ability to undergo large elastic strains without excessive perma- 
nent set. This property enables the test model to be buckled 
several times in axial compression, at virtually the same buck- 
ling load--the so called "repeatability'' phenomenon. According 
to the manufacturer, Mylar has very nearly isotropic mechanical 
properties, being made by a random calendering process. 
Mylar has been used recently in many tests of which in partic- 
ular, Seide, Weingarten and Morgan specifically investigated 
how Mylar cylinders could be made and tested. They determined 
that the material has a proportional limit of about 6,000 psi, 
a Young's Modulus in the neighborhood of 700,000 psi and a 
Poisson's ratio of approximately 0.3 [10]. The nominal vari- 
ation of the sheet thickness was found to be well within those 
limits specified by the manufacturer(+ 0.001 inches). 

The above values of Young's Modulus and Poisson's ratio 


were used for all calculations performed in this investigation. 
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3. Model Parameters 

Since the objective of the test program was to test for the 
variation of the number of circumferential buckles and the post 
buckling characteristics of cylinders with R/t below R/t equal 
to 800, a wide range of geometric parameters was chosen for 
the tests. These values are shown in the Schedule of Parameters, 
Table I. 

The thickness of the cylinders was a nominal 0.0075 inches 
(7.5 mils) or 0. 010 inches (10 mils). The radius was varied 
from two to five inches and the length from 2.75 inches to 
10.0 inches. 

4, Test Models 

All test models were thin-walled circular cylindrical shells 
with a glued lap joint along the generated line. The models were 
fabricated by wrapping Mylar sheet about a wooden cylindrical 
mandrel. The Mylar sheets were accurately cut into the desired 
axial lengths and the circumferential length of the desired cyl- 
indrical model plus one inch overlap for the bonding of the lap 
foint. Bonding was accomplished by using an uncured silicone 
rubber adhesive. To insure secure bonding, the overlap was 
restrained by a wooden bar with felt stripping, which was held 
securely against the lap joint by circumferential bands until the 


adhesive had cured. The cylindrical mandrel had a flat plate 


FZ 


attached at one end to aid in the forming of the test model 
squarely. Although the thin Mylar sheet is quite flexible, the 
use of carefully laid out patterns and wooden mandrels made 
it possible to fabricate test models that were dimensionally 
accurate and relatively free of initial wrinkles or bulges. 

End clamping was accomplished by inserting the model 
into a circular cavity of an accurately machined end plate and 
filling the cavity with rubber cement. This edge condition was 
utilized for all the tests conducted in this investigation. Great 
care was taken to insure that the end plates were parallel so as 
to obtain stress uniformity when the load was transmitted to 
the test model by means of the end plates. 

It has often been observed experimentally that a small in- 
ternal pressure will raise the buckling load of a cylinder sig- 
nificantly. To obtain comparable data it is essential to maintain 
a uniform pressure. Ventilation holes were therefore cut in 
the end plates in order to equalize internal and external pres- 
sures so that no pressure buildup would occur. 

The test model fabrication procedure is pictorially dis- 
played in Figure 3. A view of the various test model sizes is 
shown in Figure 4. 

5. Testing the Model 


The test machine used in this investigation was a standard 
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Baldwin 60,000 pound Universal Test Machine. All of the 
tests conducted were performed on this machine. The test 
models were centered accurately relative with the testing ma- 
chine. The load was applied until buckling had occurred. 

In this type of investigation it is most important that the 
Wise o is evenly distributed around the circumference of the 
test model. The alignment of the test model with the test 
machine loading plate was accomplished by the uSe of a ball 
ay This technique is the most convenient method of insuring 
that the load is applied parallel to the test mmdel axis. The 
ball joint consists of two metal circular plates in each of which 
is cut a spherical dome, and into which a greased ball is in- 
serted. Initially in the testing procedure the ball joint is 
placed at the center of the end plate where the load is applied. 

The first Pee tests conducted on a given model were used 
to adjust the eon of the ball joint so as to obtain as nearly 
as possible a symmetrical buckling pattern. Proper alighment 
was then judged to have been achieved, regardless of the po- 
sition of the ball joint through which the load was applied. 

When the test model had failed, the load was held as near 
as possible at the drop-off value while the buckles it had formed 
were measured. For 60% of the models investigated the circum- 


ference of the model was not filled completely with buckles; 


14 


the buckles that formed were measured. The wavelength was 
then used to extrapolate the buckles to the complete circumference. 
The loading and unloading cycle was repeated several times 
on each test model to obtain repeatable values of the buckling load 
and pattern. After each unloading the ball joint was displaced 
away from the direction where the previous buckling pattern had 
been initiated. This procedure helped to insure proper centering. 
A minimum of five tests was completed on each test model. 
The testing machine and model setup are shown in Figure 5. 
6. Results | 
Over 300 tests were conducted on over 50 test models to 
obtain sufficient data for the analysis. Room temperature for all 
tests was approximately 72°F, plus or minus two degrees. 
There was no control on the humidity which was assumed to have 
remained constant. 
Test model failure was characterized by the snap through 
action of the shell wall into the buckled condition. Figure 6 is 
a photograph of a typical buckling pattern observed during the tests. 
Bad results were catagorized into obviously bad (poor fabrica- 
tion rejects-two test models), and suspected bad (outside the 
population of the normal pattern of results-four test models). 
Bad tests were probably due to: obvious imperfections of the 
test model, which were visually observed, or test results 


(probably due to imperfections) 
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which statistically fell outside the population of the normal pattern 
of results. Obvious imperfections were due to poor fabrication or 
wrinkling of the test model prior to buckling. 

For the majority of tests the buckling initiated locally and it 


was necessary to extrapolate the buckle wavelength around the 


Te 


~ 


circumference of the model to obtain the number of circumferen- 
tial buckles. This condition is attributed to a non-uniform stress 
pattern in the wall of the model which is caused by an off centered 
loading or a tilted “Tne eevrousiee this results in a decrease in 
the buckling load compared to the ideal experimental condition. 

In this analysis the above condition was acceptable, since as long 
as the models tested supplied results which fell into a homogeneous 
population, they were meaningful in the buckling analysis. 

Results of the test arorapapecaune conducted are found in Tables II 
and III with pertinent relations graphically displayed in Figures 
7 eee lor 

| The first series of tests, identified as the control tests, were 
conducted on test models of the same geometry in order to de- 
termine if results would fall into the same population. Test re- 
sults are found As TableII, Summary of Data for Control Tests, 

with a graphical display of the results in Figures 7 and 8. The 
results of the control tests indicate that the frequency of distri- 


bution of the buckling load and the circumferential buckle 
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number, N, is normal and that the small scatter of results is 
about a definite meaningful value. 

From Table III it is obvious that as the testing program 
progressed the results obtained were more consistent. This is 
attributed to the learning process, whereby as the testing and 
fabrication program progressed there was more experience to 
rely upon. 

As shown in Figure 9, Variation of N with R/t and L/R, 
Theory and Experimental Results Compared, the tests confirm the 
primary objective of the analysis. The number of circumferential 
buckles for the region of R/t less than 800 is definitely de- 
termined by the geometric parameters R/t and L/R and the 
value of N is approximately predicted by de Neufville's theoret- 


1/4 
; R/t 


ical expression: 


It is noted that the test results fall very close tothose theoret- 
ically predicted. 
From Figures 10 through 15 it is seen that the variation 


of L/R alone does not change o/o ° Investigation of the 


.° 
values of o/o , obtained for the ratios of L/R chosen inthe 


analysis at a given R/t ratio shows that o/o | is independent 


1 
of the parameter L/R. This result is in agreement with a 


recent analysis by another investigator [9]. It is also observed 
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from Figures 10 through 15 that as the parameter R/t is 
increased in value, the ratio of o/o | decreases. 

‘ayeont Figure 16, Variation of o/o 4 with R/t, it is shown 
that as R/t increases, the value of o/o 4 decreases. This 
result is contrary to de Neufville's conclusion that o/o was 
independent of the parameter R/t. In Figure 16 it is noted 
that the parameter R/t has its greatest effect at ratios of R/t 


below 500. For values of the ratio above 500 the ratio of 


o/o 


26) is observed to remain nearly constant. The logical 


reason for the above effect is that for larger, thinner shells 
with a higher ratio of R/t, there are likely to be more imper- 
fections. There are more imperfections found due to the larger 


surface area. 


7. Conclusions 

This investigation adds to the general study on experimental 
effort developing general nee on the buckling of thin- 
walled cylindrical shells. The polyester film, Mylar, is adapta- 
ble and definitely a significant aid in such a study. The ease in 
the fabrication of cylindrical shells out of ee sheet, and the 
repeatability phenomenon, enabling the model to be buckled sev- 
eral times at virtually the same buckling load, enables the 
investigator to obtain considerable test data with a minimum of 


experimental effort and time. 
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The investigation conducted confirms that the number of cir- 
cumferential buckles, N, is a function of both the parameters, 
L/R and R/t. When R/t is held constant, N will decrease 
with increasing values of the parameter L/R. When L/R is 
held at a constant value, N will increase with increasing values 
of the ratio of R/t. It is important to determine the value of 
N because it is an index of the buckling mode, which in turn con- 
trols the behavior and strength of the shell after buckling. From 
the investigation it was found that N accurately follows the ap- 


proximate theoretical expression: 


1/4 
7 Rit 
wane [BE] 


developed by deNeufville. 

It is also concluded, for the range of the parameters, L/R 
and R/t, considered in the investigation, that the ratio of the 
buckling stress to that of the classical buckling stress is: 

(1) independent of the parameter L/R, and (2) dependent on 
the parameter R/t, with increasing values of the ratio of R/t 


resulting in a decrease in o/G | The parameter R/t has its 


] 

greatest effect at ratios of R/t below 500. For values of the 

ratio above 500 the ratio of o/o remains nearly constaht. 
Further experimental effort on the buckling of cylindrical 


shells using the polyester film, Mylar, can definitely expand and 


contribute greatly to the general understanding of the buckling 
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mechanism. Likely areas of investigation are: (1) larger spread 
in the variation of the parameter L/R, (2) values of the R/t 
parameter greater than 1600, (3) effect of reinforcement of 

ne shell wall by the use of circumferential stiffeners, (4) effect 
of the axial overlap used for the bonding of the lap joint, used in 
this investigation and earlier investigations by Weingarten, Morgan, 
and Seide, and de Neufville, and (5) effects:of various end fixity 
conditions (deNeufville stated that, from a theoretical point of 

view, the study of the buckling mechanism has no preferable form 


of edge restraint [7].) 
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Table I 


Schedule of Parameters 


t R L et L/R 
7.5 mils 2" 1. 83" 267 918 
| 2.51 | 1.25 
4, 0" 2.0 
3" 2.751 400 918 
| 3.75" | 1.25 
Geo!" 20 ¥ 
5M 4, 58" , 667 918 
| 6. 25" | 1.25 
10. 0" 200 
10 mils 2" 1. 83" 200 918 
| 2.5" | 1. 25 
4, 0" 2.0 
31 2.75" 300 918 
| 3.75" | 1,25 
6. 0" 2.0 
51 4, 58" 500 918 
| 6.25" | 1.25 
10. 0" 2.0 
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Table II 
Summary of Data for Control Tests 
R/t = 400 L/R= 1.58 
= 7.5 mils R= 3.0 inches 


Run Buckling load 


paoce! No. (lbs. ) 


> 


ae GC) 


= Oo wm FP W NO EF DB oO FF W ND 


N DPD wo & W WNW 
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74 
Ue 
70 
66 
66 
68 
(a0 
72 
(GZ 
69 
1 
69 
tz 
¢2 
0 
fez 
14 
74 
eZ 


9-0 
oan 
Cag 
Guo 
Ona 
920 
7. 
9.4 
9.4 
950 
90 
8.8 
oa 
oe 
§.6 
See) 
8.9 
ve (2 
vB 


Model 


Run 


- om -&_ W WH 


oO ONO Fe W LH 


Table II (Cont'd. ) 
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Buckling load 
(lbs. ) 


66 
68 
66 
62 
64 
72 
7 
75 
78 
78 
78 


9.4 
a. 9 
8.8 
2. 
2.) 
9.0 
as, 
77 
9.4 
>. 
se 


Table Ill 


Summary of Data for Test Models 
2 


. reject 
t Model Average No.of Average 

R/t L/R (mils) No. O /o, tests N 
400 . 916 1s 7 . 404 6 WOES) 
| 8 oe 6 1036 
9 . 403 5 1 Om 
25 ] . 448 y 9.6 
Z . 458 10 9% 
3 vols 8 vee) 
3A .417 7 92 
2.0 13, 588 7 8.4 
| 14 .407 7 820 
Ws) . 362 6 Sa 
300 . 916 NO) 10 . 448 8 Orc 
| 1] poor fabrication reject 
2 .378 6 9.4 
W743) “ .457 6 Seno 
| 5 ps). )9) 5 8.22 
6 . 454 7 8.8 
z2.0 16 374 9 8.1 
| 17 441 7 S60 
| 18 553 8 7.9 
267 5 DIN: C25 26 Sa 5 ie 
| 27° 653 6: 9.1 
28 . 465 a, vies 
| eras: Ze "otS 10 8.4 


Ze 


Bable 4H dCont'de«) 


ee se 
23 598 6 8.5 
| 24 473 6 8.4 
25 526 7 8.4 
2.0 1S (eas Ee fran 
| 20 563 7 7.5 
21 605 8 7.4 
200 .916 10.0 36 498 7 8.4 
| a an ae; 
38 505 7 8.5 
is 32 534 6 7.9 
33 493 7 7.9 
34 661 6 7.8 
35 622 6 7.8 
2.0 29 peee fabrication reject 
| | 30 589 8 6.9 
31 486 6 6.9 
667 946 7.5 43 417 5 11.2 
44 368 5 11.1 
1.25 41 354 6 10.5 
42 re 5 10.5 
2.0 39 417 5 9.7 
40 410 6 9.5 


fet | 


Table III (Cont'd. ) 





Bit WEIR mits) Noe ala teste NE 
500 916 O70 50 385 b 10.4 
51 2206 5 10.6 
1.25 47 351 6 o 7 
48 . 432 6 9.8 
a0 . 368 5 9.6 
2.0 | 45 . 384 6 8.8 
46 . 382 6 8.8 
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Fliigge, 1932 (theoretical) o/o,)] = 0.6Et/R 


Wilson and Newmark, Lundquist, Donnell, 1934 (Ghd re- i 
1932, Donnell, 1934 (experimen- tical) imperfection theory 
tal) Wo.) 0.4 and decreases as 

increases. 


Ballerstedt and Wagner, 1936 
Kanemitsu and Nojima, 1939 


(experimental) o/o,, also de- 


creases as L/R increases 
for small L/R. 







von Karman and Tsien, 
1941 (theoretical) demon- 
strate that a non-linear 
post-buckling curve exists 
and that low equilibrim 

' states are possible. 








Tsien, 1942 (theoretical) 
energy balance and hump. 







Batdorf, Schildcrout and Stein, Michielsen, 1948 (theore- 
1947 (curve-fitting) L/R effect tical) minimum o/o v 

f cl 
attributed to edge effects. 0.4. 


Donnell and Wan, 1950 (the- 
oretical) more refined im- 
erfection theory. 












Hoff, 1962 (theoretical) 
o/o <a 0.3Et/R fora 


cylinder with free edges. 


Kempner, 1954 (theoret- 


ical) minimum o/o , w 
a cl 







Babcock and Sechler, Tennyson, 
1963 (experimental) more per- 
fect shells, o/0 ~ 0.8. 













Almroth, 1963 (theoreti- 
cal) minimum o/o.) ~ 
Q.A.22 


Hoff, 1964 (theoretical) 
o/o., = 0.3Et/R for ellip- 


tical case N=2. 






de Neufville, 1965 (experimental) distinguishes between elastic 
and plastic buckling of cylinders; postbuckling stress depends on 
N and thus on L/R and R/t, and minimum tends to zero; buck- 
ling stress independent of L/R and R/t alone, but imperfections 
are demonstrated to increase for the longer, thinner shells. 





Figure l. 
Chronology of investigation of buckling stress 
of thin cylinders [10]. 
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Wilson and Newmark, Lund- 
quist, 1933. Donnell, 1934 
(experimental) N Ome 


von Karman and Tsien, 1941 
(theoretical) N depends on 
the shape of the buckle. For 
square buckles, n = OQ. 91VR/t 
N decreases as end shorten- 
ing increases. 


Donnell, 1950 (theoretical) 


“N= 0.91 V R/t. Observed 
N is one half of N predicted. 


Thielemann, 1962 (experimen- 
tal) observed that N decreases 
as end shortening increases. 
Also reports secondary buck- 
ling phenomenon. 


de Neufville, 1965 (experi- Forsberg, 1964 (theore- 


mental) N increases as R/t tical) Vibration mode 
increases, decreases as analysis: N increases 
L/Rincreases: — 1/4 as R/t increases, de- 
Niceemer Bayt creases as- L/R in- 
| L/R creases. 


Obtains postbuckling stress- 
strain path for several N. 





Figure 2. 
Chronclogy of the investigation of the 
buckling mode of thin cylinders [10]. 
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Ficure 3, Test model fabrication procedure 
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Figure 10. w/t) versus L/R for R/t equal to 200 
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Figure 11. TAZ} versus L/R for R/t equal to 267 
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Figure 12. T/T, versus L/R for R/t equal to 300 
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Figure 13, Tigss versus L/R for R/t equal to 400 
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Figure 14. Ca versus L/R for R/t equal to 500 
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Figure 15, oko versus L/R for R/t equal to 667 
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